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In this paper we study the self-propulsion of a triangular micro-robot (or A-robot) 
which consists of three spheres connected by three rods; the rods' lengths are chang- 
ing independently and periodically. Using the asymptotic procedure containing the two- 
timing method and distinguished limit arguments, we obtain analytic expressions for 
sclf-propulsion velocity the angular velocity. Our calculations show that a A-robot ro- 
tates with constant angular velocity around its centroid, while the centroid moves in a 
circle. The important special case of zero angular velocity represents rectilinear transla- 
tional sclf-propulsion with constant velocity. 



1. Introduction 



The studies of micro-robots represent a flourishing modern research topic which strives 
to create a fundamental base for modern applications in medicine and technology, see e.g. 
Purcell (1977)| |Becker et.al. (2003)[ [Najafi fc Golestanian (2004)1 preyfus et.al. (2005) 



Earl et.al. (2007)|| Chang et.al. (2007)1 |Alouges etaL (2008)[[Gilbert at.al. (2010)HGolestanian & Ajdari (2008) 
Golestanian fc Ajdari (2009)[|Alexander et.al. (2009)[|Leoni et.al. (2009)[|Lauga (2011) 



Romanczuk et.al. (2012) The simplicity of micro-robot's geometry represents the ma- 
jor advantage in contrast to extreme complexity of self-swimming microorganisms, e.g. 
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Pedley fc Kessler (1987) Vladimirov et.al. (2004) Pedley (2009) Polin et.al. (2009) This 



advantage allows us to describe the motion of micro-robots in greater depth. 

In this paper, we study the self-propulsion of a triangular three-sphere micro-robot (we 
call it A-robot), see the figure. We formulate the problem for a A-robot with arbitrary 
sides, but the final expressions are presented only for the case when all mean sides are 
equal. Our calculations show that the A-robot rotates with constant angular velocity Q 
around its centroid, while the centroid itself moves along a circle with the same angular 
velocity J7. The important special case of f2 = represents rectilinear translational self- 
propulsion with constant velocities. To obtain these results we employ the two-timing 
method and distinguished limit arguments, which lead to a simple and rigorous analytical 
procedure. Our approach is technically different from all previous methods employed in 
the studies of micro-robots (except Vladimirov (2012b)| Vladimirov (2012c) ). The possi- 
bility to describe explicitly any motion of a A-robot shows the strength of our method. 



The used version of the two-timing method has been developed in Vladimirov (2005) 
Vladimirov (2008 )[| Vladimirov (2012a) 
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2. Formulation of problem 

We consider the A-robot consisting of three rigid spheres of equal radii R connected by 
three rods of different lengths k, I, and m. In plane cartesian coordinates (x, y) the centers 
of the spheres are = (x^\ ), v = 1,2,3 and the triangle (with vertices in centres 
of spheres) is composed of vectors k = x^ 1 ^ — x^ 3 \ I = x^ — x^\ and m = x^ — x^ 2 \ 
see the figure. We describe the centers of the spheres as 

x (1) = X + a, x {2) = X + b, x {3) = X + c 

where X = (X, Y) is a centroid radius- vector such that vectors a, b, c, unit vectors a, 
b, c, and angles 9, ip, \ are introduced as 



a + b + c = 0; k = a — c, I = b a, m = c b; 



a = aa, b = bb, c = cc 



(2.1) 
(2.2) 



cos6> » 



i sm u j 



I cos(6» + V) ' 



sin(0 + ip) 



' cos(6» + ip + x) 



c = 



sin 



4> + x) 



Well-known relations between the sides and medians of a triangle yield 



k 2 = 2(a 2 + c 2 ) - b 2 , I 2 = 2(a 2 + 6 2 )-c 2 , m 2 = 2(b 2 + c 2 ) - a 2 , (2.3) 



see Johnson (2007) The A-robot moves due to the prescribed oscillations of the sides 



k = k = k + efc(r), 1 = 1 = 1 + d(r), m = rh = m + em(r); r = (jt (2.4) 

where u> = const e = const; k, I, m are mean lengths, and k, I, rh are 27r-periodic functions 
of t with zero average values. The spheres experience external friction forces F^ u ' = 
(F[ v , F% ), while the rods are so thin in comparison with i?, that their interaction with 
a fluid is negligible. The considered problem contains three characteristic lengths: the 
length of rods L, the radius of spheres R, and the amplitude of rod's oscillations a. The 
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characteristic time-scale is T and the characteristic force is F. We have chosen 

3L=k + J + m, a = eL, T = l/u, F = 6ttt]RL/T (2.5) 

where 77 is viscosity of a fluid. Two independent small parameters are 

s=a/L, <5 = 3i?/(4L) (2.6) 

The dimensionless variables (marked with asterisks) are chosen as 

x = Lx\ t = Tt*, Fi = F F* ; (2.7) 

Below we use only dimensionless variables, however all asterisks are omitted. 

The form of a triangle is completely prescribed by three lengths I, k, m or, alternatively, 
by three lengths a, b, c (|2.2[) . (|2.3j) . We choose generalised coordinates as 

<? = (91,92,53,94,95,36) = (X,Y,0,a,b,c) (2.8) 

The choice a, b, c as generalised coordinates drastically simplifies further calculations. 
The motion of the A-robot is described by the Lagrangian function C = C(q, q t ), which 
includes constraints (|2.4[) with Lagrangian multipliers a, /3, 7 

q t )=K + a(a - a) + P(b -b)+ 7 (c - c) (2.9) 

where subscript t stands for d/dt, given functions a, b, c correspond to I, I, rh in (|2.4j) . (|2.3[) , 
and K, is kinetic energy of a robot. Lagrangian multipliers represent additional unknown 
functions of time. The Lagrange equations are 

dtdqnt oq n dq n 
where Q = (Qi,Q2,Q3) is the generalized external force, exerted by a fluid on the A- 
robot. As one can see, we use latin subscripts (i,k = 1,2) for cartesian components of 
vectors and tensors, subscript n for generalised coordinates, and subscripts (or super- 
scripts) H,v = 1,2,3 to identify the spheres. The fluid flow past the A-robot is described 
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by the Stokes equations, where all inertial effects are neglected. Since masses of spheres 

and rods arc negligible, then K = 0. Hence, (f231 . (|2.10j) , (j23j) . ([2~2"j) produce the system 

of equations 

+F {2) +F {3} = (2.11) 
a x F (1) + b x F (2) + c x F (3) = (2.12) 
a-F (1) = -aa, b ■ F {2) = -(3b, c-F (3) =- 7 c (2.13) 
a = a, 6 = 6, c = c (2-14) 

where the subscript _L stands for a perpendicular to (x, y)-plane component of a vector. 
One can see that for chosen generalised coordinates the Lagrangian equations appear as 
zero total force (|2.11|) and zero total torque (|2.12[) . The explicit expressions for are 

F n ~ _ z m + 6 n^\ = J2 s^x^ (2.15) 

Each force represents the first approximation for the Stokes friction force exerted 
on a sphere moving in the flow generated by other two spheres. To construct (|2.15[) we 



use a classical explicit formula for fluid velocity past a moving sphere, see Lamb (1932) 



Landau fc Lifshitz (1959)[|Moffatt (1996)| The equations (|2TTT|) - ([2~T5)) represent a system 
of nine equations for nine unknown functions of time: X, Y, 9, a, 6, c, a, fi, and 7. It is 
well known from analytical mechanics that Lagrangian multipliers represent the reactions 
of constraints; hence the equations (|2.13p can be kept out of consideration if we are not 
interested in the forces exerted by rods. 

Since the explicit form of system f|2.1L() - f)2.15|) is rather cumbersome, we restrict our- 
selves to the special case of an equilateral mean triangle k = 7 = m = 1. Using (|2.1[) we 
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can rewrite the correspondent equations (|2.11|) . (|2.12|) as 



3X t = S 



§ (12 >(2X t - c t ) + S (13) (2X t - bt) + § (23) (2X t - a t )] (2.16) 



[axa t + bxb t + cx c t ] ± =i|ax n (1) + b x n (2) + c x II (3) J ^ (2.17) 

3. Two-timing method and asymptotic procedure 

3.1. Functions and notations 
The following dimensionless notations and definitions are in use: 

(i) s and r denote slow time and fast time; subscripts r and s stand for related partial 
derivatives. 

(ii) A dimensionless function, say G = G(s,t), belongs to class I if G = 0(1) and 
all partial s-, and r-derivatives of G (required for our consideration) are also O(l). In 
this paper all functions belong to class X, while all small parameters appear as explicit 
multipliers. 

(iii) We consider only periodic in t functions {G E V : G(s, r) = G(s, r + 27r)}, where 
s-dependence is not specified. Hence, all considered below functions belong to V {\L. 

(iv) For arbitrary G € V the averaging operation is 



1 



T0 + 27T 



(G> = — y G(s,T)dT = G(s), Vr (3.1) 

(v) The tilde-function (or purely oscillating function) represents a special case of V- 
function with zero average (G ) — 0. The bar-function (or mean-function) G = G(s) does 
not depend on r. A unique decomposition G = G + G is valid. 

3.2. Asymptotic procedure and successive approximations 

The introduction of fast time variable r and slow time variable s represents a crucial 
step in our asymptotic procedure. We choose r = t and s = e 2 t. This choice can be 



justified by the same distinguished limit arguments as in Vladimirov (2012a) Here we 
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present this choice without proof, however the most important part of this proof (that 

this choice leads to a valid asymptotic procedure) is exposed and exploited below. We 

use the chain rule 

d/dt = d/dT + s 2 d/ds (3.2) 

and then we accept (temporarily) that r and s represent two independent variables. In 
further consideration we consider series in small parameter e and restrict our attention 
to the terms 0(e 2 ). Simultaneously we keep at most linear terms in S. In our theory S 
appears not separately but as products such as e 2 8. We do not specify the dependence 
of unknown functions on 6; such dependence reveals itself naturally during calculations. 
Unknown functions are taken as regular series 

X(r, s) = X {t, s) + eX x {T, s) + e 2 X 2 (r, s) + . . . (3.3) 

and similar series for 9, a, b, c, a, (3, and 7. We take 

X = and 6 = 0, while X ^ and 9 ^ 

which express the basic property of our solutions: long distances of self-swimming and 
large angles of rotations are caused by small oscillations. We also have accepted ao = 
~b~o = Co = 1/ v3 for an equilateral triangle. After the application of p.2[) to Q3.3P we have 

X t =eX lT + e 2 (x 2T + X 0s ) + O{e 2 ) (3.4) 

and similar expression for 8. In calculations below all the bar functions belong to the zero 
approximation and all the tilde-functions belong to the first approximation, therefore we 
omit the related subscripts. 

The successive approximations of equations (|2.16[) . (|2.17p yield: 
Terms O(e ) give identities = 0. 



Terms O(e 1 S ) lead to 
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X T = 0, 39 T + 2ip T + Xr = 



(3.5) 



where the second equation follows from 



axa T + bxb T + cxc =0 



(3.6) 



and definitions (|2.2[) . The integration of (|3.5[) in the class of periodic functions yields 



X = 0, 36» + 2ip + x = 



(3.7) 



From the law of cosines one can derive that xp = a + b — 2c and x — b + c — 2a, then (|3.7[) 



»1\'(>S 



= c-6. 



(3.8) 



Terms O(eS): These terms do not vanish and can be easily calculated. However, (as one 
can see below) they do not participate in the leading terms of average motion that appear 
in the order 0(e 2 ). 
Terms 0(e 2 ): Eqns. (|2T6|) . (|2~T7)) give 



[a x a s + b x b s + c x c s ] 



a x a T + b x b T + c x c T 



(3.9) 
(3.10) 



From these equations one can already see that the leading terms Xo s = 0(8) and 9q s 
0(1). In further calculations all tilde-functions, say, 

" — (Sik + 3lilk)lr + hhr + hkr, etc. 



ikr 



(3.11) 



should be expressed in terms of a, b, c with the use of (|2.1j) . (|2.2[) . Then straightforward 
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but rather cumbersome transformations of p.9j) yield 



3X s /6 = (ab T )c + (ca T )b + (bc T )a 

a/3 ~ 1 ~ ~ 

= —^-(ab T — ca T )a + — (2bc T — ab T — ca T )a 

4 ~_ 

S = —={ac T + bc T + c T b) 
V3 

where unit vectors a , 6 , and c correspond to ^-derivatives of a, 6, and c, say, 



(3.12) 



(3.13) 



a = 



V 



sin ( 



cos( 



b = — (— a + \/32'), b = — (— \/32 — a'), e£c. 



One can check that expressions ()3.12j) . (|3.13[) possess a relabeling invariance (a,b,c) i-> 
(b,c, a), eic. The final expressions are 



X s = -(Aa + Ba'), 9 S = D: 
y 

A = — ((/ + k)rh T ), B = -(2kl r - lrh r - rhk T ) 



D = - 



3a/3 



(I m T + kl T + m k T ) 



(3.14) 
(3.15) 
(3.16) 



with constants A, A', and il. It is also useful to exclude m by introducing the angle <p 
between the sides k and I. The law of cosines gives m = (fc + / + \/3^/2)/2; its substitution 
to (i3~T5l) . (i3T6)) yields 

A=l{(l + k)<p T ), Bs^ + ^I-*)^), D = -i((r-fc)^ T )(3.17) 

Both equations (|3. 14=|) can be immediately integrated as 

8 S 
X-X = —(Asin9 + Bcos8), Y -Yo= — (-Acosd + Bs'md) (3.18) 

6 = Ds + 6> n 



with constants of integration Xq, Yq, and 9q. Finally, we should recall that s is slow time 
variable, hence physical translational velocity and angular velocity are 



(5p 2 

V = X t = (Aa + Ba), n = 6 t =e 2 D 

9 



(3.19) 
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For dimensional velocities V and f2 one has to add multipliers Lw and ui correspondingly 
(see (|2.5p ). It is noticeable that V and have different amplitudes: V = 0(Se 2 ), but 
SI = 0(e 2 ). It corresponds to an important physical fact: translational self-propulsion of a 
A-robot takes place due to hydrodynamic interactions between spheres, while rotational 
self-propulsion does exist without taking such interactions into account. Eqns. (|3.18[) show 
that the A-robot rotates with constant angular velocity J! around its own centroid and 
simultaneously moves with the same angular velocity f2 in a circular orbit of radius 
SyJ A 2 + B 2 /9il and centered at the point (Xq, Yq), which can be defined from the initial 
data. 



4. Examples of rectilinear self-propulsion 

The condition for motion p,14|) to be purely translational is f2 = 0; it could be really 
useful for practical applications. The related absolute value of translational velocity is 



|V| = 5e 2 \/A 2 + B 2 /9 (4.1) 

Hence, the equation 

(lfh T + kJ T + fhk T ) = or ((I - k)(p T ) = (4.2) 

provides options for purely translational (rectilinear) motion. The simplest cases corre- 
spond to isosceles oscillations, when k = I ^ fh or k = fh ^ I or I = fh ^ k. For example, 
k = I rh gives 

il = 0, A = V3{km r ) =3(T^ T )/4, B = 



which coincides with the translational self-swimming of a F-robot by Vladimirov (2012c) 



However, in Vladimirov (2012c) the mean angle between the arms of a y-robot is arbi- 



trary, while everywhere in this paper it is 7r/3. Hence, the results of Vladimirov (2012c) 



are complementary to the study done in this paper. It is also worth to emphasise that 
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variables and the method of calculation in Vladimirov (2012c) are different, hence the 
obtaining of the same result strongly suggests that the calculations of this paper do not 
contain mistakes. 

Another case of rectilinear motion is <p = 0, which corresponds to a fixed angle between 
the sides k and I, while k ^ and 1^0. From ()3.17[) one can see that in this case A = 0, 
D = 0, and B = 3(lk T )/2. It gives velocity V = Se 2 (lk T )a' /6, which shows that the 
A-robot is not rotating but propagating in direction m. This example can be seen as 



generalization of the micro-swimmer by Golestanian & Ajdari (2008) for a fixed angle 



7r/3 between the arms. Similarly, the examples of pure rotation (without translation) can 
be build. 



5. Discussion 

(i) We have derived that for any oscillations of three sides of the considered A-robot 
the averaged paths of each sphere represent circles with centers that move in circles 
(|3.18[) . That is, the paths of each sphere are epicycloids. Following other examples in 



fluid dynamics (see Bennett (2006) ) one can call these paths Ptolemaic solutions. 

(ii) The magnitude 0(e 2 S) of translational velocity in terms of small parameters 
(|3.19|) . (|4~Tj) is the same as the result by |Golestanian fe Ajdari (2008)[|Golestanian fe Ajdari (2009)[ 



Vladimirov (2012b) for linear micro-robots. At the same time, our choice of slow time 



s = e 2 t (|3.2j) agrees with classical studies of self-propulsion for low Reynolds numbers, 



see 



Taylor (1951)] |Blake (1971)] |Childress (1981)| as well as the geometric studies of 



Shapere fe Wilczek (1989) 



(iii) The magnitude 0(e 2 ) for angular velocity ()3.19j) represents a novel result. It 
corresponds to an important physical fact: rotational self-propulsion docs exist without 
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taking hydrodynamic interactions between spheres into account; it caused by standard 
Stokes' drag force in an infinite fluid and the reactions of constraints. 

(iv) We have built an asymptotic procedure with two small parameters: e — > and 
S — > 0. Such a setting usually requires the consideration of different asymptotic paths on 
the plane (s,S) when, say S = 5(e). In our case we can avoid it, since 5 does not appear 
separately, but only in combinations such as e 2 5. 

(v) The mathematical justification of the presented results can be performed by the es- 
timation of an error in the original equation, as in |Vladimirov (2010)||Vkdimirov (2011) 



It is also possible to derive the higher approximations of V and fi, as it has been done 
by |Vladimirov (2010)| |Vladimirov (201 1)| for different cases. These approximations can 



be useful to study the cases when self-propulsion in the main order vanishes. 

(vi) In order to compare velocities of micro-robots and micro-organisms we use di- 
mensional variables, where V ~ ujLe 2 6; it shows that a A-robot can move itself with 
the speed about 10% of its own size per second (we have taken e = S = 0.2 and u> = 



30s , the value of u> can be found in Pedley & Kessler (1987) Vladimirov et.al. (2004) 



Pedley (2009)] |Polin et.al. (2009) | From these papers we also can see that this estima- 



tion of self-propulsion velocity is about 10 times lower than a similar value for natural 



micro-swimmers. 



The author is grateful to Profs. M. Bees, A. Gilbert, K.I. Ilin, and J. Pitchford for 
useful discussions. 
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